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Abstract 



We consider a parametrically forced pendulum with a vertically oscillating 
suspension point. It is well known that, as the amplitude of the vertical os- 
cillation is increased, its inverted state (corresponding to the vertically-up 
configuration) undergoes a cascade of "resurrections," i.e., it becomes stabi- 
lized after its instability, destabilize again, and so forth ad infinitum. We make 
a detailed numerical investigation of the bifurcations associated with such res- 
urrections of the inverted pendulum by varying the amplitude and frequency 
of the vertical oscillation. It is found that the inverted state stabilizes via 
alternating "reverse" subcritical pitchfork and period-doubling bifurcations, 
while it destabilizes via alternating "normal" supercritical period-doubling 
and pitchfork bifrucations. An infinite sequence of period-doubling bifurca- 
tions, leading to chaos, follows each destabilization of the inverted state. The 
critical behaviors in the period-doubling cascades are also discussed. 
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I. INTRODUCTION 



We consider a parametrically forced pendulum whose suspension point undergoes a verti- 
cal periodic oscillation. The system is described by a second-order non-autonomous ordinary 
differential equation [Q], 

le + be + ml{g - euj'^cosujt) sine = 0, (1) 

where the overdot denotes the differentiation with respect to time, / is the total moment of 
inertia, 6 is a damping coefficient, m is a mass attached to one end of a light rigid rod (its 
mass can be negligible) of length /, 6 is the angular displacement measured counterclockwise 
from the downward vertical, and e and u are the driving amplitude and frequency of the 
vertical oscillation of the suspension point, respectively. Making the normalization ut — > 27rt 
and 6 — >■ 2nx, we obtain a dimensionless normalized form of Eq. (|l]), 

X + 2'Kj3Vtx + 2-K{n^ - A cos 27rt) sin 2tcx = 0, (2) 

where ujq = '^mgl/I, /? = b/Iuo, Q = Uq/u, and A = mle/I. 

The parametrically forced pendulum has two stationary states. One is the "normal" 
state corresponding to the vertically-down configuration with x = 0, and the other one is 
the "inverted" state corresponding to the vertically- up configuration with x = |. For the 
case of the "unforced" simple pendulum (with A = 0), the normal state is obviously stable, 
while the inverted state is clearly unstable. However, as the normalized driving amplitude A 
is increased above a critical value, the inverted state becomes stable. This stabilization of the 
inverted pendulum has been discussed theoretically and demonstrated experimentally 

uni- 

Here we are interested in the bifurcations associated with stability of the inverted state. 
As in the case of the normal state [O , the linear stability of the inverted state is determined 



by a damped Mathieu equation |jT3| 



ii + 2tt pnii + ATT^{-n^ + A cos2Trt)u = 0. (3) 
It is well known that the Mathieu equation has an infinity of alternating stable and unstable 



A ranges for a given f2 Consequently, as the parameter A is increased, the inverted 

pendulum exhibits a cascade of "resurrections" (i.e., it becomes stabilized after its instabil- 
ity, destabilizes again and so forth ad infinitum) for any given Q. By varying the normalized 
driving amplitude A and the normalized natural frequency Q, we make a detailed numeri- 
cal investigation of bifurcation behaviors associated with such resurrections of the inverted 
pendulum for a fixed value of the normalized damping coefficient (3. 

This paper is organized as follows. In Sec. |T|, we discuss bifurcations associated with sta- 



bility of periodic orbits, using the Floquet theory |15]. The bifurcation behaviors associated 



with the resurrections of the inverted state are then investigated through numerical calcula- 



tiona of its Floquet (stability) multipliers in Sec. |IT|for the case (3 = 0.2. It is found that the 
stabilizations of the inverted state occur via alternating "reverse" subcritical pitchfork and 
period-doubling bifurcations, while its destabilizations take place through alternating "nor- 
mal" supercritical period-doubling and pitchfork bifrucations. After each destabilization 
of the inverted state, an infinite sequence of period-doubling bifurcations follows and leads 



to chaos. In Sec. IV, we also study the critical behaviors in the period-doubling cascades. 



Finally, a summary is given in Sec. 0. 
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II. STABILITY, BIFURCATIONS, LYAPUNOV EXPONENTS, AND WINDING 

NUMBERS 



In this section, we first discuss stability of periodic orbits in tlie Poincare map of tlie 
parametrically forced pendulum, using the Floquet theory ||15| . Bifurcations associated with 



the stability, Lyapunov exponents and winding numbers are then discussed. 

The second-order ordinary differential equation (|^) is reduced to two first-order ordinary 
differential equations: 

x = y, (4a) 
y = f{x, y, t) = -27ipny - 27r((]^ - A cos 2nt) sin 27rx. (4b) 

These equations have the (space) inversion symmetry S, because the transformation 

S -.x ^ -X, y ^ -y, t^t, (5) 

leaves Eq. (^ invariant. If an orbit z{t)[= {x{y),y{t))] is invariant under S, it is called a 
symmetric orbit. Otherwise, it is called an asymmetric orbit and has its "conjugate" orbit 
Sz{t). 

The surface of section for the parametrically forced pendulum is the Poincare (time- 
1) map. Hence the Poincare maps of an initial point zq [= {xQ,yo)] can be computed by 
sampling the orbit points Zm at the discrete time t = m {m = 1,2,3,...). We call the 
transformation Zm — > Zm+i the Poincare map and write Zm+i = P{zm)- 

The linear stability of a g-periodic orbit of P such that P'^{zq) = zq is determined from the 
linearized-map matrix DP"^ of P'^ at an orbit point zq. Here P'^ means the g-times iterated 
map. Using the Floquet theory |]T5|, the matrix DP'^ can be obtained by integrating the 



linearized differential equations for small perturbations as follows. 

Let z*(t) = z*(t+q) be a solution lying on the closed orbit corresponding to the g-periodic 
orbit. In order to determine the stability of the closed orbit, we consider an infinitesimal 
perturbation 6z[= {u,v)] to the closed orbit. Linearizing Eq. (^ about the closed orbit, we 
obtain 




1 

fx{X , t) fy 



(6) 



Here and fy denote the partial derivatives of f{x,y,t) in Eq. (|) with respect to the 
variables x and y, respectively. They are given by 

(x, t) = -Att^ {n^ - A COS 2nt) cos 27rx, fy = -27r/3fi. (7) 

Note that J is a 2 x 2 g-periodic matrix. Let W{t) = {w^{t),w'^{t)) be a fundamental solution 
matrix with W{0) = I. Here w^{t) and w'^{t) are two independent solutions expressed in 
column vector forms, and I is the 2x2 unit matrix. Then a general solution of the g-periodic 
system has the following form 



3 



Substitution of Eq. (||) into Eq. (|^) leads to an initial-value problem to determine W{t) 

W{t) = J{t)W{t), W{0) = I. (9) 

It is clear from Eq. (|^) that W{q) is just the linearized-map matrix DP'^{zq). Hence the 
matrix DP'^ can be obtained through integration of Eq. (P) over the period q. 
The characteristic equation of the linearized-map matrix M(= DP'^) is 

A2-trMA + detM = 0, (10) 

where trM and detM denote the trace and determinant of M, respectively. The eigenvalues, 
Ai and A2, of M are called the Floquet (stability) multipliers. As shown in [jT6[, detM is 
calculated from a formula 

detM = ei'o'f-^'^*. (11) 



Substituting the trace of J (i.e., trJ = —2t[[3VL) into Eq. (pJ]), we obtain an exact analytic 
result 

det M = (12) 

(Note that detM is a constant, independently of the orbits.) Accordingly, the pair of Floquet 
multipliers of a periodic orbit with period q lies either on the circle of radius e~^^^'^ or on 
the real axis in the complex plane. The periodic orbit is stable only when both Floquet 
multipliers lie inside the unit circle. We first note that they never cross the unit circle, 
except at the real axis and hence Hopf bifurcations do not occur. Consequently, a stable 
periodic orbit can lose its stability when a Floquet multiplier decreases (increases) through 
— 1 (1) on the real axis; conversely, an unstable periodic orbit can gain its stability when a 
Floquet multiplier increases (decreases) through —1 (1) on the real axis. 

When a Floquet multiplier A decreases through —1, the stable periodic orbit loses its 
stability via period-doubling bifurcation. On the other hand, when a Floquet multiplier 
A increases through 1, it becomes unstable via pitchfork bifurcation. For each case of the 
period-doubling and pitchfork bifurcations, two types of supercritical and subcritical bifur- 
cations occur. For the supercritical case of the period-doubling and pitchfork bifurcations, 
the stable periodic orbit loses its stability and gives rise to the birth of a new stable period- 
doubled orbit and a pair of new stable orbits with the same period, respectively. On the 
other hand, for the subcritical case of the period-doubling and pitchfork bifurcations, the 
stable periodic orbit becomes unstable by absorbing an unstable period-doubled orbit and a 
pair of unstable orbits with the same period, respectively. Hereafter, all these bifurcations, 
associated with instability of a stable periodic orbit, will be called the "normal" bifurcations. 
We also note that reverse processes of the normal bifurcations can occur for the case of un- 
stable orbits. That is, when a Floquet multiplier of an unstable orbit increases (decreases) 
through —1 (1), it becomes stabilized via "reverse" period-doubling (pitchfork) bifurcations. 
For example, for the reverse subcritical period-doubling and pitchfork bifurcations, the un- 
stable orbit gains its stability by emitting an unstable period-doubled orbit and a pair of 



unstable orbits with the same period, respectively. For more details, refer to Ref. |17|. 

We now discuss the Lyapunov exponent and the winding number of an orbit in the 
Poincare map P. Expressing the linearized equations (|^) for the displacements in terms of 
the polar coordinates u = r cos (p and v = r sin 0, we have 
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r = r[(l + /^) sine/) COS + /y sin^ </)], (/) = — sin^ (/) + (/i^ cos</) + sin <;/>) cos<^. (13) 

The motions of the displacements (r, 0) contain all the information about the nearby orbits. 
Hence we first obtain the Poincare maps of an initial displacement (tq, 0o) by sampling the 
displacements (r^, 0m) at the discrete time t = m {m = 1,2,3,...). Then the average 
exponential rate of growth of the radius r, 

a = lim —In—, (14) 

gives the largest Lyapunov exponent a, characterizing the average exponential rate of diver- 
gence of the nearby orbits. If a is positive, then the orbit is called a chaotic orbit; otherwise, 
it is called a regular orbit. On the other hand, the average rate of increase of the angle 
(normalized by the factor 27r), 

uj = lim (15) 

m^oo 27rm 

gives the winding number w, characterizing the average rotation number of the nearby orbits 
during the time 1 (i.e., one iteration of P). For more details on the Lyapunov exponent and 
the winding number, refer to Ref. |jTI 



III. BIFURCATIONS OF THE INVERTED STATE AND TRANSITIONS TO 

CHAOS 

In this section, by varying the two parameters A and Q, we study bifurcations associated 
with stability of the inverted state for a damped case of /3 = 0.2. It is well known from 
the theory of the Mathieu equation |jl^ that there exist an infinity of disconnected stability 



regions in the fl — A plane (i.e., an infinity of alternating stable and unstable A ranges 
exist for any given Q). Consequently, as A is increased, the inverted state undergoes a 
cascade of resurrections (i.e., it stabihzes after its instability, destabihzes again, and so forth 
ad infinitum) for any given Q. We make a detailed numerical investigation of bifurcation 
behaviors associated with such resurrections of the inverted state. 

As explained in Sec. ^ the linear stability of a periodic orbit with period q in the 
Poincare map P is determined from the linearized-map matrix M(= DP'') of P'^. The 
matrix M can be obtained through numerical integration of Eq. (^) over the period q, and 
then its eigenvalues give the Floquet multipliers of the periodic orbit. In such a way, we 
determine the stability regions of the inverted state in the Q — A plane through numerical 
calculations of its Floquet multipliers A's. The first three stability regions, denoted by Sn 
{n = 1,2,3), are shown in Fig. |l|. Each Sn is bounded by its lower stabilization curve, 
denoted by L„, and by its upper destabilization curve, denoted by f/„. As the order n is 
increased, the stability region Sn becomes smaller. 

We investigate bifurcation behaviors associated with the stabilizations and destabiliza- 
tions of the inverted state at the stability boundary curves in detail. It is found that they 
depend on whether the order n of the stability region Sn is odd or even. At the stabilization 
curves L„ of odd (even) n, the unstable inverted state becomes stable via reverse pitchfork 
(period-doubling) bifurcations. However, the situation becomes reverse for the case of the 
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destabilization curves Un- That is, the stabihzed inverted state loses its stabihty through nor- 
mal period-doubling (pitchfork) bifurcations when the destabilization curves of odd (even) 
n are crossed. These period-doubling and pitchfork bifurcation curves are denoted by the 
solid and dashed curves in Fig. |l|, respectively. There are two types of supercritical and 
subcritical bifurcations for the case of the period-doubling and pitchfork bifurcations, as 
explained in Sec. |Il|. All the stabilizations occur via the subcritical bifurcations, while all 
the destabilizations take place through the supercritical bifurcations. Consequently, with 
increasing A the inverted state stabilizes via alternating reverse subcritical pitchfork and 
period-doubling bifurcations, while it destabilizes through alternating normal supercritical 
period-doubling and pitchfork bifurcations. After each destabilization of the inverted state, 
an infinite sequence of period-doubling bifurcations, leading to chaos, also follows and ends 
at a finite accumulation point. We obtain such accumulation points for several values of Q 
{fl = 0.001, 0.05, 0.1, 0.15, 0, 2). They are denoted by solid circles in Fig. |I] and seem to lie 
on smooth critical lines. When crossing these critical lines, period-doubling transitions to 
chaos occur. 

The lower and upper bifurcation curves, L„ and Un, in Fig. |1| are also labelled by the 
winding numbers u of the inverted state as Ln{uj) and Un{uj), respectively. We obtain the 
winding number u of the inverted state through numerical integration of the linearized 
equation (13) over the period 1. It is known that for the pitchfork bifurcations, the winding 
numbers of the inverted state become integers, while they are odd multiples of 1/2 for the 



period-doubling bifurcations [|T8[. Note that the winding number u of the inverted state 
increases with respect to the order n of the bifurcation curves. 

We now present the concrete examples of bifurcations associated with the resurrections 
of the inverted state for the case Q = 0.1. The bifurcation diagrams and the phase- flow 
and Poincare-map plots are also given for clear presentation of the bifurcations. We first 
investigate the bifurcation behaviors associated with the first resurrection of the inverted 
state with increasing A. For the unforced case of A = 0, the inverted state is clearly 
unstable. However, when the lower stabilization curve Li of the first stability region 5*1 
is crossed at the first stabilization point As{l) = 0.142 066 . . ., the unstable inverted state 
becomes stabilized via reverse subcritical pitchfork bifurcation. The bifurcation diagram 
near the first resurrection of the inverted state is shown in Fig. @(a). Through the reverse 
subcritical pitchfork bifurcation, a conjugate pair of unstable asymmetric orbits with period 
1 appears, and their phase portraits for A = 0.15 are shown in Fig. H(b). However, when 
the upper destabilization curve Ui of S'l is crossed at the first destabilization point Ad{l) = 
0.471 156 . . ., the stabilized inverted state loses its stability via normal supercritical period- 
doubling bifurcation. Consequently, a stable period-doubled symmetric orbit appears and its 
phase portrait for A = 0.5 is also shown in Fig. |^(c). Note that the winding number u of the 
inverted state increases from to 1/2, as A is changed from ^4^(1) to ^4^(1). We also study the 
subsequent bifurcations with increasing A further. Unlike the case of the inverted state, the 
symmetric 2-periodic orbit becomes unstable by a symmetry-breaking pitchfork bifurcation, 
which leads to the birth of a conjugate pair of stable asymmetric orbits with period 2. 
Then each 2-periodic orbit with broken symmetry undergoes an infinite sequence of period- 
doubling bifurcations, ending at its accumulation point A^ (= 0.575 154 . . .). Consequently, 
period-doubling transition to chaos occurs when the parameter A increases through Al. The 
critical scaling behaviors of period doublings near the critical point Al are the same as those 
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for the ID maps, as will be seen in Sec. |I^. 

With further increase of A, we also study the bifurcations associated with the second 
resurrection of the inverted state. Since the order of 5*2 is even, the types of bifurcations 
associated with the stabilization and destabilization become different from those for the case 
of the first resurrection. When the lower stabilization curve L2 of S2 is crossed at the second 
stabilization point ^4^(2) (= 3.779 771 . . .), a reverse subcritical period-doubling bifurcation 
occurs, which is in contrast to the case of the first resurrection. Consequently, the unstable 
inverted state stabilizes with the birth of a new unstable symmetric orbit with period 2. 
Figure |^(a) shows the bifurcation diagram in the vicinity of this second resurrection. The 
phase portrait of the newly-born unstable symmetric 2-periodic orbit for A = 3.783 is also 
shown in Fig. ^(b). However, when the upper destabilization curve U2 of S2 is crossed 
at the second destabilization point Ad{2) = 3.811973 . . ., a normal supercritical pitchfork 
bifurcation occurs, which is also in contrast to the case of the first destabilization. As a 
result, the stable inverted state destabilizes with the birth of a conjugate pair of stable 
asymmetric orbits of period 1. The phase portraits of the newly-born stable asymmetric 
orbits with period 1 for A = 3.815 are shown in Fig. ^(c). Note that the winding number lj 
also increases from 1/2 to 1, as A is varied from ^4^(2) to Ad{2). As A is further increased, 
a second infinite sequence of period-doubling bifurcations, leading to chaos, also follows and 
ends at its accumulation point A2 (= 3.829 784 . . .). The critical scaling behaviors of period 
doublings near A = A2 are the same as those near the first accumulation point A*. 

Finally, we investigate the bifurcations associated with the third resurrection of the 
inverted state. The types of bifurcations associated with the stabilization and destabilization 
become the same as those for the case of the first resurrection, because the order of 5*3 is odd. 
As shown in Fig. Q^a), the unstable inverted state stabilizes with the birth of a conjugate 
pair of unstable asymmetric orbits with period 1 via reverse subcritical pitchfork bifurcation 
at the third stabilization point ^^(3) (= 10.671579 . . .) on the lower stabilization curve 
L3. As A is increased, the stable inverted state also destabilizes with the birth of a stable 
period-doubled symmetric orbit through a normal supercritical period-doubling bifurcation 
at the third destabilization point ^4^(3) (= 10.673 818 . . .) on the upper destabilization 
curve U3. The phase portraits of the newly-born orbits for the cases of the stabilization and 
destabilization are shown in Figs, ^(b) and (c) for A = 10.6718 and 10.6741, respectively. 
We also note that, as A is changed from As{3) to Ad{3), the winding number uj increases 
from 1 to 3/2. Since the winding number of the inverted state increases, its nearby orbits 
have an increasing number of loops [e.g., compare Figs. |^(b) and (c) with Figs. |^(b) and 
(c)]. The subsequent bifurcation behaviors are also the same as those for the above first 
case. Consequently, a third infinite sequence of period-doubling bifurcations, leading to 
chaos, follows and ends at its accumulation point ^3 (= 10.675 090 . . .). The critical scaling 
behaviors of period doublings are also the same as those near A^. 

IV. CRITICAL BEHAVIORS IN THE PERIOD-DOUBLING CASCADES 

In this section, we first investigate the winding-number sequence of the period-doubling 
cascade and find that the winding numbers at the period-doubling bifurcation points con- 
stitute an alternating sequence converging to a limit value, as in other oscillaotrs [0. The 
orbital scaling behavior and the power spectra of the periodic orbits born via period-doubling 
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bifurcations as well as the parameter scaling behaviors are then investigated. These critical 



scaling behaviors for all cases studied are found to be the same as those of the ID maps [|19 

As an example, we consider the case Q = 0.1. The first three accumulation points 
y4*'s {i = 1,2,3) of the period-doubling bifurcations are shown in Fig. Only the critical 
behaviors at the first accumulation point Al are given below, because the critical behaviors 
at all the three accumulation points are the same. For this first case, we follow the periodic 
orbits of period 2^ up to level k = 8. As explained in Sec. |T|, the stabilized inverted state 
loses its stability for A = Ad{l) through a normal supercritical period-doubling bifurcation, 
giving rise to the birth of a stable period-doubled symmetric orbit. However, this symmetric 
2-periodic orbit becomes unstable via symmetry-breaking pitchfork bifurcation, which results 
in the birth of a conjugate pair of asymmetric orbits with period 2. Then each 2-periodic 
orbit with broken symmetry undergoes an infinite sequence of period-doubling bifurcations, 
ending at its accumulation point Al. Therefore, period-doubling transition to chaos takes 
place when the parameter A increases through A^, 

Figure |^(a) shows the bifurcation diagram of the first period-doubling cascade. [For the 
sake of convenience, only one asymmetric orbit of period 2 is shown in Fig. ^(a).] The largest 
Lyapunov exponent a and the winding number u in the period-doubling cascade are also 
given in Figs. |§(b) and (c), respectively. The largest Lyapunov exponent has a constant value 
(= —irjSQ) when the Floquet multipliers A's lie on the circle of radius e~'"^^'^ (in the complex 
plane), while it changes smoothly when A's lie on the real axis. The value of a becomes 
zero at each period-doubling bifurcation point. Unlike the case of a, the winding number lo 
takes a constant rational value when A's lie on the real axis, and hence it can change only 
when A's lie on the circle of radius e~'^^^^. Consequently, oj[A) becomes a step-like function. 
Rational steps appear near the period-doubling bifurcation points, as shown in Fig. ^(c). 
The sequence of the winding-number steps is given by tUfc = |(1 — -^-^f-) {k = 0, 1, 2, . . .); 
the first three values (corresponding to k = 1,2 and 3) are given in Fig. §(c). Note that 
the winding numbers constitute an alternating sequence converging to its limit value uj^ 
(= 1/6). Consequently, the quasi-periodic attractor at the first accumulation point A* has 
the winding number uj^d- study the winding-number sequences in the second and 

third period-doubling cascades. They are given by cufc = |(2 + ^f") and a;^ = |(7 — ^^^^ ), 
respectively. As in the first case, these alternating sequences also converge to their limit 
values uj^ (= 2/3) and cu^ (= 7/6). Hence, the quasiperiodic attractors at the second and 
third accumulation points A2 and A^ have their winding numbers a;^^ and o;^-' , respectively. 
Note that the winding numbers cu^ of the quasiperiodic attractors at the accumulation 
points A* increase with i. 

Table | gives the A values at which the period-doubling bifurcations occur; at A^, a 
Floquet multiplier of an asymmetric orbit with period 2'^ becomes —1. The sequence of A^ 
converges geometrically to its limit value Al with an asymptotic ratio 6: 

4 = - 5. (16) 

The sequence of 6^ is also listed in Table I. Note that its limit value 6 (~ 4.7) agrees well with 
that (= 4.669 ■ ■ ■) for the ID maps 0. We also obtain the value of A* (= 0.575 154 232 . . .) 
by super converging the sequence of {Ak} pOl) . 

As in the ID maps, we are also interested in the orbital scaling behavior near the most 
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rarified region. Hence we first locate the most rarified region by choosing an orbit point 
^(fc) (^a^C^)^ ^C*^))] which has the largest distance from its nearest orbit point (z^''"') for 
A = Ak. The two sequences {x^'^^} and {y^^^}-, listed in Table |T|, converge geometrically to 
their limit values x* and y* with the ID asymptotic ratio a (= —2.502 ■ ■ ■), respectively: 

^ _ ^(fc) "^'^ ^ y{k+l) _ y{k) ~^ (-'-''') 

The values of x* (= 0.350 686 . . .) and y* (= 0.014 623 . . .) are also obtained by supercon- 
verging the sequences of x*^^-* and y^''\ respectively. 

We finally study the power spectra of the 2^-periodic orbits at the period-doubling 
bifurcation points A^. Consider the orbit of level k whose period is g = 2*^, {z!;li> = 
(^m^2/m'*)' ^ — 0) 1) ■ ■ ■ ) 9 ~ !}• Then its Fourier component of this 2'^-periodic orbit is 
given by 

z^'\u;,) = -j:zi:;^e-^-^-, (18) 

y m=0 

where uj = 2nj /q, and j = 0, 1, . . . , g — 1. The power spectrum P^^\ujj) of level k defined 
by 

P^'^\u,) = \z^'\u,)\\ (19) 

has discrete peaks at a; = Uj. In the power spectrum of the next {k + 1) level, new peaks 
of the {k + l)th generation appear at odd harmonics of the fundamental frequency, ujj = 
27r(2j + l)/2('=+^) (j = 0, ... ,2^ - 1). To classify the contributions of successive period- 
doubling bifurcations in the power spectrum of level k, we write 

P^'^ = PooSico) + E E ' - ^O), (20) 

1=1 j=0 

where Pip is the height of the jth peak of the Ith generation appearing at = ujij (= 
27r(2j + l)/2'). As an example, we consider the power spectrum P^^''{u) of level 8 shown in 
Fig. p. The average height of the peaks of the Ith generation is given by 

'^^^HO = ^'e"<- (21) 

It is of interest whether or not the sequence of the ratios of the successive average heights 

2/3('=) (/)=(/)('=) (/)/(/)('=) (/ + 1), (22) 

converges. The ratios are listed in Table |T|. They seem to approach a limit value, 2(3 ~ 21, 
which also agrees well with that (= 20.96 ■ ■ ■) for the ID maps [^ . 
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V. SUMMARY 



We make a detailed investigation of bifurcations associated with resurrections of the 
inverted state through numerical calculations of its Floquet multiphers. It is found that 
its stabilizations occur via alternating reverse subcritical pitchfork and period-doubling bi- 
furcations, while its destabihzations take place through alternating normal supercritical 
period-doubhng and pitchfork bifurcations. An infinite sequence of period-doubUng bifur- 
cations, leading to chaos, also follows each destabilization of the inverted state. The orbital 
and parameter scaling behaviors near the accumulation points A* of the period-doubling 
cascades are also found to be the same as those of the ID maps, although the winding 
numbers uj^ of the quasi-periodic attractors at the accumulation points increase with i. 
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TABLES 

TABLE L Asymptotically geometric convergence of the parameter sequence {^/t} 



k Ak 6k_ 

1 0.573 847 671035 

2 0.574 992 231118 8.26 

3 0.575130 862 947 7.25 

4 0.575149 971647 5.68 

5 0.575153 335 892 4.78 

6 0.575154039 822 4.65 

7 0.575154191110 4.69 

8 0.575154 223 350 



TABLE II. Asymptotically geometric convergence of the orbital sequences and 



1 


0.356 951938 




0.014184 290 




2 


0.349 165 197 


-4.016 


0.014680 027 


-6.189 


3 


0.351 104216 


-3.547 


0.014 599 925 


-2.709 


4 


0.350 557 529 


-3.078 


0.014 629 491 


-3.543 


5 


0.350 735119 


-2.608 


0.014 621 145 


-2.409 


6 


0.350 667034 


-2.509 


0.014 624 610 


-2.615 


7 


0.350 694 173 


-2.502 


0.014623 285 


-2.448 


8 


0.350 683 324 




0.014623 826 





TABLE III. Sequence 2(3^''\l) [= <j)^^\l) / + 1)] of the ratios of the successive average 



heights of the peaks in the power spectra 



k 


4 


5 


I 


6 


7 


6 


34.7 


24.6 








7 


34.0 


24.5 




21.7 




8 


33.9 


24.0 




21.7 


21.5 
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FIGURES 



FIG. 1. Stability diagram of the inverted state in the parametrically forced pendulum. The first 
three stability regions of the inverted state, denoted by Sn {n = 1, 2, 3), are shown in (a) and (b). 
For each Sn, a period-doubling bifurcation occurs on the solid boundary curve, while a pitchfork 
bifurcation takes place on the dashed boundary curve. The lower and upper bifurcation curves, 
Ln and C/„, are also labelled by the winding numbers u of the inverted state as Ln{uj) and Un{oj), 
respectively. The accumulation points of an infinite sequence of period-doubling bifurcations, 
denoted by solid circles, seem to lie on the smooth critical lines. For other details, see the text. 

FIG. 2. (a) Bifurcation diagram (plot of x versus A) in the vicinity of the first resurrection of 
the inverted state. Here the solid and dashed lines denote stable and unstable orbits, respectively, 
and q denotes the period of an orbit, (b) Phase portraits for A = 0.15. The stabilized inverted 
state is denoted by the solid circle. On the other hand, the phase flows of a conjuate pair of 
unstable orbits with period 1 born via reverse subcritical pitchfork bifurcation are denoted by 
dashed curves, and their Poincare maps are represented by the crosses, (c) Phase portraits for 
A = 0.5. The destabilized inverted state is denoted by the cross. On the other hand, the phase 
flow of a stable orbit with period 2 born via normal supercritical period-doubling bifurcation is 
denoted by a solid curve, and its Poincare maps are represented by the solid circles. 

FIG. 3. (a) Bifurcation diagram near the second resurrection of the inverted state. The solid 
and dashed lines and q represent the same as those as in Fig. 1. (b) Phase portraits for A = 3.783. 

The stabilized inverted state is denoted by the solid circle. On the other hand, the phase flow of 
an unstable orbit with period 2 born via reverse subcritical period-doubling bifurcation is denoted 
by a dashed curve, and its Poincare maps are represented by the crosses, (c) Phase portraits for 
A = 3.815. The destabilized inverted state is designated by the cross. On the other hand, the phase 
flows of a conjugate pair of stable asymmetric orbits with period 1 born via normal supercritical 
pitchfork bifurcation are denoted by solid curves, and their Poincare maps are represented by the 
solid circles. 

FIG. 4. (a) Bifurcation diagram near the third resurrection of the inverted state. The solid 
and dashed lines and q denote the same as those in Fig. 1. Note that the bifurcation behaviors 
associated with the stabilization and destabilization of the inverted state are the same as in Fig. 1. 
The phase portraits of the orbits associated with the stabilization and destabilization are shown 
in (b) and (c) for A = 10.6718 and 10.6741, respectively. For other details see the text. 

FIG. 5. (a) Bifurcation diagram [plot of x versus — ln(^J — A)] of the first period-doubling 
cascade, (b) plot of the largest Lyapunov exponent a versus — ln(A^ — ^4), and (c) plot of the 
winding number u versus — ln(AJ — A) for O = 0.1. Here {q, s or a) denotes the stable A range of 
the symmetric or asymmetric orbit with period q. 

FIG. 6. Power spectrum P(8) (a;) of level 8 for A = As (= 0.575 154 223 350) . 
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